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Abstract 

Group theoretical technique for construction of two-electron wavefunctions with zero total mo- 
mentum in solids based on the Mackey theorem on symmetrized squares is developed. These states 
correspond to singlet and triplet Cooper pairs. The nodal structure of these states is investigated 
group theoretically and compared with experimental data for unconventional superconductors. It 
is shown that when the Mackey theorem is applied twice, the main four-electron states in solids 
can be constructed. Possible connection of these states with experimental data is discussed. 
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I. INTRODUCTION 



The antisymmetrical many-electron wavefunctions constitute the basis for theory of 
strongly correlated electronic systems in atoms molecules and sohds The mathematical 
basis for such a technique in the spherically symmetrical case is well developed (see e.g. 
[1,2]). On the other hand the applications of permutational symmetry methods to solid 
state wavefunctions are rather scant. The method of two-electron wavefunctions construc- 
tion in solids is mainly due to Bradley and Davies [3] who made use of Mackey theorem on 
symmetrized squares of induced representations [4] . Theory of relatively recently discovered 
unconventional superconductivity (i.e. high Tc [5] and heavy fermion [6]) is not developed 
yet, and symmetry considerations are very useful for understanding of electron structure of 
these new materials. Non-totally symmetrical order parameter (or the wavefunction of a 
Cooper pair) with line and point nodes on Fermi surface is the mail feature of unconventional 
superconductors [7]. The induced representation method"^ was used for construction of an- 
tisymmetrical wavefunctions in solids of different symmetry. [8-12]. This group theoretical 
approach revealed a new symmetrical feature of two-electron functions in solids : their nodal 
structure. Comparison of group theoretical nodal structure for Dq^ group with experimental 
data for unconventional superconductor U Pt^ resulted the symmetry for superconduct- 
ing order parameter (SOP) [11]. The space group approach makes possible construction of 
two-electron wavefunction by projection operator technique [12]. But it is clear from the 

— # — * 

shape of these functions that the commonly accepted interaction for a Cooper pair in k — k 
manifold cannot give different energies for different irreducible representations (IRs) . Hence 
it follows that the interactions between other wave vector are important in solids. Group 
theoretically it means that one should construct the four electron wavefunctions with zero 
total momentum i.e. Cooper quartets. In the present work the space-group approach to 
two electron wavefunction is briefly reviewed and the theory is generalized to four-electron 
wavefunctions. The classification of four electron states connected by space inversion and 
mirror reflection is given making use Young tables 
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II. THO-ELECTRON WAVEFUNCTION IN SOLIDS (COOPER PAIRS) 



Prom a unitary IR tk of a subgroup H one can construct a unitary representation of the 
whole group G [13]. The structure of this unitary representation (induced representation) 
depends on the left coset decomposition of the whole group with respect to its subgroup: 

G^J2s,H (1) 

i 

where i — l....n and n — \G\/\H\ 

The induced representation is defined by the following formula : 



{tn T G){g)i^j^ = t^{Si ^gsj)i,j{si ^gsj, H) (2) 
1, if sj^gsj e H 

where:5(Si gSj.H) = { 

0, if Si gsj ^ H 

Following [13] we use an up directed arrow for the notation of induction. The indexes i 
and j in formula denote the block columns and rows of the induced representation matrix 
and correspond to the single coset decomposition (1). The indexes /i and v number the rows 
and columns of the "small" IR tk- 

In the case of crystal symmetry the induced representation (2) is irreducible representa- 
tion of a space group, provided the group if is a wave vector k group ( little group) and tk 
is its unitary IR (small IR). The action of left coset representatives Sj on the wave vector k 
results all prongs of its star |^|. In the case strong spin-orbit coupling the IRs tk in formula 
(3) are replaced by double- valued small IRs Pk ■ 

According to the Pauli exclusion principle the total two-electron wavcfunction is antisym- 
metric with respect to permutation of electronic coordinates. Hence in a weak spin-orbit 
coupling [L — S scheme) the symmetrized Kronecker square of the spatial part of the wave- 
function is combined with antisymmetrized Kronecker square of its spin part (singlet pair) , 
and the antisymmetrized Kronecker square of the spatial part of the wavefunction is com- 
bined with the symmetrized Kronecker square of its spin part (triplet pair). In a strong 
spin-orbit coupling case [j — j scheme) the wavefunction belongs to the antisymmetrized 
Kronecker square or double-valued IR of the space group. According to Anderson [14] the 
Cooper pair wavefunction is invariant with respect to lattice translations Hence it follows 
the consideration is limited by the centre of a Brillouin zone for two-electron states. 
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The structure of the Kronecker square of an induced IR may be envisaged by the double 
coset decomposition of G relative to H which is written (3) as: 

G = J2Hd^H (3) 

(7 

The sum runs over all distinct double cosets a. 

Corresponding wave vector k is defined by the following formula 

k+dak = ka+ba (4) 
The intersection of wave vector groups in the left hand side is written as: 

M^^Hn d^Hd-^ (5) 

For each double coset we consider a representation of subgroup M^- defined by the formula: 

= t,{m) X t,{d-^md,) (6) 

Where m G Mg- 

For self-inverse double coset, i.e. : 

Hd^H = Hd-^H (7) 

there are two extensions of Pa into the subgroup: 

M« = M« + (8) 

where a — d^hi — h2d^ and /ii, /i2 £ H. 

These extensions corresponding to symmetrized and antisymmetrized parts of Kronecker 
square are defined in terms of their characters as follows: 

x(Pc^(am)) = +x(tn{,amam)) (9) 
x{P~{am)) = -xit^iamam)) (10) 

where m e M^- 
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The symmetrized and antisymmetrized parts of the Kronecker square of induced repre- 
sentation are written by two following formulae respectively (the Mackey theorem [3] on 
Kronecker squares): 

[t« T G X T G] = [q, X q,] T G + ^P+ ^G + Y^P^ T G (11) 
{t« T G X T G} = {t, X T C + T G + T G (12) 

a fj 

The first items on the right-hand sides of (11) and (12) correspond to the double coset 
defined by the identity element, a corresponds to self- inverse double cosets and (3 to non- 

self-invcrse double cosets for which HdpH ^ Hd^^H. In the case of a strong spin-orbit 
coupling case the possible symmetries of two-electron states are obtained by substituting of 
double- valued IRs into formula (12). 

— * 

If the one-electron wave vector k belongs to a general point inside a Brillouin zone the 
two-electron wave vector, defined by formula (4), equals zero for the self- inverse double coset 
defined by the space inversion. The extended intersection group Ma defined by formula (8) 
is the group Gj consisting of two elements : E and /. For single-valued IRs we immediately 
obtain that Pa equals to IR Ag of group Gj and that Pa equals to IR of group Gj. 

These representations are induced into the central extension of the space group (point 
group) . The induced representation can be easily decomposed making use of the Probenius 
reciprocity theorem: the number of appearance of the IR of the whole group in the 
decomposition of the induced representation pk ^ G equals to the number of appearance 
of IR Pfc-iii the decomposition of F^^when it subduccd to the subgroup. Making use of 
Frobenius theorem we obtain two conclusions. Firstly, in agreement with Anderson (14) , 

— * 

we obtain that for k a general point of a Brillouin zone all even IRs are possible for singlet 
pairs and all odd IRs are possible for triplet pairs. Secondly, the number of appearance of 
each IR equals to its dimension. Hence it follows that for one-dimensional IR the result is 
unique, but for two-dimensional IRs there are two non-equivalent basis functions and one 
can take any linear combinations. From this point of view experimentally observed double 
superconducting transition in UPts (15) may be connected with two non-equivalent states 
corresponding to the same two-dimensional IR . To obtain total wavefunction of a Cooper 
pair in a weak spin-orbit coupling {L — S scheme) one should multiply the spatial part of 
the wavefunction by spin singlet function 5"° for singlet pair and by spin triplet function 
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for triplet pair. 

In strong spin orbital coupling case the representation Pa equals to IR of group Cj 
and even IRs are missing. To obtain all possible pair symmetries, the time reversal 6 should 
be considered. In the absence of magnetic fields the total symmetry of a crystal with Fedorov 
group G is described by the Shubnikov II (grey) magnetic group: 

M^G + eC (13) 

where ^ is a time-reversal operation 

The time-reversal symmetry results in additional degeneration for one- dimensional small 
double- valued IRs i.e. at general points and at the planes of symmetry in a Brillouin zone. To 

obtain all possible two-electron states one should use induced corepresentations D{pk) T G 
[13] in formula (12). 

— * 

For k a general point in a Brillouin zone the decomposition of corepresentation P~ (see 
formula (6)) contains representations Ag and 3^4^ of the group Cj. The IR Ag corresponds 
to singlet pair and "iAu correspond to three components of triplet pair. 

The superconducting state is usually more ordered than normal state, i.e. the transition 
to it is accompanied by the symmetry reduction[16] . One possible way is the time-rcvcrsal 
symmetry violation, i.e. transition from the direct product 6 xG to ordinary Fedorov group 
G or to one the Shubnikov group 9 x {G — H) + H [13]. Total number of different cases 
of construction of Shubnikov groups is quite large. In order to envisage general trends 
we consider the simplified case of time reversal symmetry violation and its influence on 
the nodal structure of superconducting order parameter. Ferromagnetic fluctuations can 
be approximated as time-reversal symmetry violation. In this case the one-electron states 
belong to double valued IR of the space group. For k a general point of a Brillouin zone 
we obtain for two-electron states two IRs A^ of group : one for spin up states and one for 
spin-down state. Two remaining IRs correspond to antiferromagnetic pairs: Ag for singlet 
pair and Au for triplet pair . Hence we obtain following formula for the character of the 
possible Cooper pair representation (reducible) in the antiferromagnetic state 

Xaitif Xnormal "^Xferro (-^^) 

Where Xnormai Xferro obtained respectively by substitution of double valued small 
corepresentation and double-valued I into the formula antisymmetrized Kronecker square of 
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induced representation. It should be noted that formula (14) is valid in the case of one- 
dimensional double-valued small IR for one-electron states, i.e. at general points and at the 
planes of symmetry in one-electron Brillouin zone. 

Making use of formula(14) for a general point in a Brillouin zone we obtain that the 
Kroneckcr product decomposition for antiferromagnctic state contains IRs Ag and of the 
group Cj. Hence it follows that in antiferromagnctic state even and odd Cooper pairs are 
possible. This general result agrees with the experimental data which show both even and 
odd symmetry for antiferromagnctic heavy-fermion superconductors CeCu2Si2 and UPt^ 
respectively [7]. 

The space-group approach to the wavefunction of a Cooper pair makes it possible to in- 
vestigate the nodal structure of SOP as follows. One should consider the distinct directions 
and planes of symmetry in a one-electron Brillouin zone and calculate the antisymmetrized 
Kronecker squares with zero total momenta of double- valued IRs or of double- valued corep- 
resentations. The absence of any IR in this square indicates a node of the SOP of this 
symmetry. There are two types of nodes. The intersection of the direction of nodes of any 
IR with the Fermi surface results in the point node. The intersection of the plane of nodes 
with the Fermi surface results in the line of nodes. 

Possible IRs for all states at the plain of symmetry (group C2h) are presented in Table 
1. In normal state all odd IRs are present and one even IR Bg of the group C2/1 is absent in 
the decomposition. Hence, it follows that in this case, only nodes of even order parameter 
on the planes of symmetry are required by the space-group symmetry and no limitations on 
odd IRs exist. This statement is in agreement with the Blount [17]' theorem according to 
which, it is vanishingly improbable for "triplet" superconductors to have curves of vanishing 
gap on the Fermi surface. If the time-reversal symmetry is violated, the antisymmetrized 
square of the double- valued IR equals to IR A^ of the group C2h- Thus, in the ferromagnetic 
state, only odd IRs are possible for the SOP on the planes of symmetry. The lack of the 
second odd IR S„ signifies that some of the odd IRs of the point group are forbidden on the 
planes. The intersection of the plane with the Fermi surface results in the hue node of the 
odd SOP. 

Going over to the antiferromagnctic state we see in Table 1 that one even IR Ay (the 
same as in normal state) and one odd IR appear in the decomposition. Hence it follows 
that the symmetry requirements for line nodes of even IRs are the same as in normal state. 
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but the lines of node of odd IRs differ from that in ferromagnetic state. 

Hence it follows that the theory is in agreement with the above mentioned experimental 
on the SOP symmetry in unconventional superconductors, i.e. antiferromagnetic supercon- 
ductors may be either even (singlet) and odd (triplet) with lines of nodes. 

Another reason for violation of Blount theorem is due to crystal symmetry lower then 
Oh-Th and T^. In the case of Oh symmetry spin function belong to three dimensional IR 
Tig. Following relation is vahd for the Kronecker product of the induced IR F of the whole 
group : 

F X (P„- T G) = (F i M„ X T G (15) 

Hence we obtain that for Oh symmetry all odd IRs of the subgroup C2h are possible for 
triplet pair and Blount theorem is fulfilled. For and D^h symmetry the Mg—l and -1 (or 
Sx and Sy) components belong to IR Eg and Ms=0 (or S^) belong to IR A2g. It is natural 
to expect that due to interactions of spins with crystal field the energies of spin states Eg 
and A2g are different and only one of them corresponds to superconducting state. Since not 
all IRs of group C2h are present in the decomposition for both cases and the lines of nodes 
appear. Thus another symmetry reason for violation of Blount theorem is the lower crystal 
symmetry. 



Table 1. 

The decomposition of representations P~ and Pa 
for the planes of symmetry (group 6*2/1) 



state 


character 


decomposition 




E 


o-h 


/ 


C2 


IRs 


normal 


4 





-2 


2 


Ag + 2Ay_ + Bu 


ferromagnetic 


1 


-1 


-1 


1 




antiferromagnetic 


2 


2 








Ag + B^ 


Pi 


1 


1 


1 


1 




Pa 


1 


1 


-1 


-1 


Bu 


Pa XTig i C^^ 


3 


-1 


-3 


1 




Pa XEgl C2/1 


2 


-2 


-2 


2 


2Au 


Pa xA2gi C2/^ 


1 


1 


-1 


-1 


Bu 
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The wavefunction of a Cooper pair for the whole group may be obtained by applying 

— * — * 

projection operators technique to the singlet and triplet wavefunction in A; — A; manifold. We 
will construct these functions taking symmetry groups and D^^ of high temperature 
superconductors. Let us denote ki the wave vector chosen in the representation domain of a 
Brillouin zone. Making use of Kovalev's [18] notation /i25 for the space inversion the spatial 
parts of singlet and triplet functions are written as: 

+ (16) 

^ - (17) 

Where the superscript of denotes the number of electronic coordinate and subscript of ip 
the prong of the /c- vector star. 

Acting by /i2 (180° rotations around the axis X) on the functions (16) and (17) we obtain 
two other basis functions: 

$2 = ^2^26 + ^M (18) 
=V'2V'26-V'26V'2 (19) 

Note that in Kovalev's^ notations for Oh group multiplication of pure rotation element 
by / corresponds to adding 24 to the element number. To construct full basis for group 
we need also functions $3*^*^ and which are obtained from $1^*^ by the action of 180° 
rotations around the axes Y and Z respectively. In addition, for D4^h group the elements 
hi3, hi6 (180° rotations around the axes (110) and (110)) and hu, /ii5 (90° and 270° 
counterclockwise rotation around Z axis) are required. 

These functions span the space two-electron wavefunctions under the action of all point 
group operations. Since the space inversion is already included in the basis functions, their 
total number equals to the half of number of point group operations. The action of pure 

— * 

rotations on the initial vector ki result in a star, whose number of prongs is half of the 
number of prongs in the the wave vector star. The action of the space inversion on the basis 
vector corresponding to any prong doesn't change a vector but introduces multiplier —1 for 
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the triplet case. Making use of standard projection operator technique and functions ^{'l^ 
we easily obtain the basis functions for Cooper pairs belonging to all IRs of group. The 
results are presented in Table 2. 

Before going to the projection for D^h group it is useful to remind the following corre- 
spondence of IRs in the subduction D^h i '■ o-nd Bi — >■ Ai, A2 and B2 — > Bi, 
E — > B2 + B^. The basis functions for one-dimensional IRs of D^h group are immediately 
obtained by projection operator technique. Since each of IRs -E'g(u) appear twice in the Kro- 
necker product decomposition, there are two independent basis sets labeled by additional 
quantum numbers . Bearing in mind the above reduction scheme , we begin with basis sets 
corresponding to IRs B2 and B^ of group we obtain the remainder results of Table 2, 
denoted by superscripts a and (3 respectively. The results for triplet pairs for group 
pairs are not presented in the Table 2. To obtain wavefunctions of triplet pairs one should 
replace subscripts to w in the first column and all superscripts s to i in the second column 
without changing of the signs. 

— * 

For k a general point in a Brillouin zone all IRs are possible for Cooper pair. But when 
the /c-vector approaches any mirror plane, the mirror reflection image of k also approaches 
the /c-vector. Total number of states decreases and lines of nodes are eventual. There are 
two possibilities. If two-electron function is unchanged under the action of the reflection, 
the function under consideration is nonvanishing on the mirror plane. On the other hand, 
if the function changes its sign, two mirror counterparts are cancelling on the plane. This 
corresponds to the line of nodes. Note, that the space inversion changes the sign of the 
spatial part of the triplet function. Making use of the above rules we can easily obtain 
nodal structure of basis functions of one- dimensional IRs of groups D2h and D^h presented 
in Table 2. 

TABLE 2. Spatial parts of Cooper pair wavefunctions for point groups D2h and 1)4^ 
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D2fe 


D4h 


IR 


n^iiT'inp' fnnrtion 


IR 


DMirinp fnnrtion 






^Ig 


(J) J _l_ (J) ^ _l_ d) ^ 1 (J) « 1 (T) « 1 (T) 

^lT^"±'2T^^3T^^4T^^i3T^^14T^^i5T^^16 




^1 ^2 ^3~^4 


^2g 


^1 ^2 ^3~^4 ^13~^14~^15 ^16 


'->2g 




Bi„ 


^1~^2~^3~^4 ^13 ^14 ^15 ^16 


^dg 


^l'^2 ^3 ^4 


r>o„ 


^1 ^2 ^3^ 4^^13 14 ^15^ 16 




$*+$2 + $3 + $4 




'^'l3-'^'?5 + *14-*16 


Blu 








B211 








Bsu 









Two dimensional IRs appear twice for k a general point in a Brillouin zone. In this case 
a direct analysis of nodal structure of basis functions of Table 2 is required. The analysis 
shows that basis functions of IR vanish in the planes(lOO) and (001) and that of IR 
vanish in the planes (010) and (001). Linear combinations of these basis functions Eg ^ E^ 
vanish in planes (110) and (110) respectively and both vanish in plane (001). Hence it 
follows that only lines of nodes in basal plane follow unambiguously from the symmetry. It 
should be noted, that point group approach also results different nodal structure of different 
two-dimensional IRs [19] . 

The analysis of broad set of experimental data on the of high-Tc superconductors [20] led 
the most of the authors to the conclusion of singlet pairing and Ag SOP symmetry in these 
compounds. Angular resolved photoelectron spectra of high-Tg superconductors reveal a 
strong trough in the diagonal of xy plane indicating dx2-y2 -pairing with line of nodes. On 
the other hand some experiments reveal also totally symmetric s -pairing without nodes. 
In many cases an interplay between these two types of pairing [20] both belonging to Ag 
IR exists. It is seen from Table 2, that Ag pairing function, obtained group theoretically 
is noddles and that other IRs have nodes in the coordinate planes only. Hence it follows 
that nodal structure of high-Tc superconductors is more complex then that which follows 
from the symmetry only. To explain this one can consider two wave vectors ka and A;^ 
symmetrical with respect to diagonal of the deformed square. Note, that the orthorombicity 
[{b — a)/{b + a)] of YBCuO is about 2 % only [20]. Two types of basis functions of Cooper 
pairs belonging to Ag IR and are easily obtained from the Table 2 by introducing 
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additional subscripts a and (3. One can suppose that due to the interaction two self- vectors 
are linear combinations of these basis states : 

*f = Ca ($f,« + + $f,« + $t) + {<^l0 + + $f,^ + <^l^) (20) 

^1 = {^ia + ^la + ^la + ^t) " {^1^ + + ^1,^ + ^1^) (21) 

Both combinations belong to IR Ag of group D2h- First one corresponds to the noddles 
s-pairing and the second to the rf2:2_j^2— pairing with line of nodes in the diagonal of xy- 
plane. In the limit of zero orthorombic distortion the symmetry group is D^h and — C/j 
, subscripts a and /3 are dropped and the sums in second brackets in right hand sides of 
(20) and (21) are written as ^fg + + ^fg + $fg. In this case combination (20) belongs 
to IR Aig and combination (21) belongs to IR of the symmetry group D^h- Hence it 
follows that the nodal structure of SOP in high- Tc superconductors is defined by hidden 
symmetry D^h. 

III. FOUR-ELECTRON WAVEFUNCTIONS 

It is clear from the shape of two-electron wavefunctions presented in Table 2 that if the 

— * — * 

interaction within the k — k manifold only, is anticipated, the energies of all IRs of the 
same parity are the same. This statement contradicts with the experimental fact that the 
superconductivity in unconventional superconductors is defined by only one IR of crystal 
point group. Hence it follows that we must suppose additional residual interaction, say 
between electrons connected by mirror reflection operation cr^. In this case one should 
consider four electron state connected by space inversion, 180° rotation and mirror reflection. 
We consider the the four -electron states which may be constructed on the basis of the 

— * — * 

singlet and triplet pairs in A; — A; manifold, whose spatial part correspond to Young tables 
[2] and [1^] respectively. In general case it is done by a plethism operation [1]. In our 
case the induced structure of IRs makes possible using of the Mackey theorem as follows. 
The spatial parts of all singlet pairs are given by second term in the right hand side of 
(11). All four electron states constructed from singlet pairs are contained in the Kronecker 
product (Pa '\ G) X {Pa t G). Due the induced structure of these representations one 
can use the Mackey theorem ones more, taking any element except the space inversion as 
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Four-electron states based on singiet poirs 
orbit spin 



12 13 



3 4 2 4 




k4 



a) Singiet pair, triplet 
with respect to C2 



b) Singlet pair, triplet 
witin respect to 



FIG. 1: 



a coset representative. In this second step the symmetrized square of singlet spatial part 
is excluded, since it results the Young table [4] whose conjugate spin part [1^] is forbidden. 
Taking 180° rotation around axis perpendicular to the plane and mirror reflection as a coset 
representatives in double coset decomposition (3) we obtain all four electron states shown 
in figures 1 a) and 1 b) respectively. 

In this figure the parity with respect to k is indicated by the sign and the spin direction 
by an arrow. Note that the Young tables corresponding to both pictures are the same, but 
indexes 3 and 4 are interchanged in pictures. All these states are odd with respect to mirror 
reflection, but no limitations are imposed on the spin orientation of mirror counterparts. 

Similar procedure can be applied to triplet pairs. In this case in the second step of 
application of Mackey theorem one can make both symmetrization and antisymmetrization. 
But in the case of totally antisymmetric spatial part one obtains a four electron state with 
totally symmetrical spin part, i.e. with parallel spins. Such a particle with total spin equals 
to 2 is not appropriate candidate for superconductivity and is excluded from the present 
consideration. 

The four electron states constructed from triplet pairs by symmetrization with respect 
to double coset defined by 180° rotation and mirror reflection are shown in figures 2 a) and 
2 b) respectively. These states belong to IRs and of the group ah can have different 
parity with respect to mirror reflection. The electrons connected by mirror reflection have 
different spin directions in both cases. 
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Four electron states based on triplet pairs 



orbit 



spin 




0) Triplet pair, odd 
with respect to 




b) Triplet pair, even 
witti respect to a^^ 



FIG. 2: 



The projection operator technique for two-electron states developed in the previous sec- 
tion results all even IRs for singlet states and all odd IRs for triplet case. When the four- 
electron states are constructed, the singlet states, even with respect to 180° rotation are 
excluded. It should be noted that additional degeneration due to time reversal symme- 
try should be taken into account for two-electron states and we can apply or not apply 
time-reversal in addition to mirror reflection. It is seen from Figure 2 that in "triplet" four- 
electron case we have only one possibility for spin orientation connected by mirror reflection. 
Hence it follows that the total number of IRs possible for four-electron quasiparticles is less 
then that in the case of two-electron quasiparticle . 



IV. CONCLUSION 



The method for construction of two- and four- electron wavefunction in solids based on 
Mackey theorem on symmetrized squares of induced representations is developed. The nodal 
structure of two-electron states is investigated group theoretically. It is shown that lines of 
nodes always exist in a singlet case and appear in a triplet case if several symmetry operations 
are violated. It is shown that the main four-electron states (related to both singlet and triplet 
pairs in k — k manifold) are described by Young table [2, 2]. Two advantages for introducing 
of four-electron states into the theory of superconductivity are emphasized. Firstly they 
correspond to interactions beyond k — k manifold and are able to explain different energies 
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for different singlet and triplet IRs. Secondly total number or permitted symmetries for 
SOP is reduced. 
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